Two-orbital Kondo effect in quantum dot coupled to ferromagnetic leads 
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We study the Kondo effect of a two-orbital vertical quantum dot (QD) coupled to two ferromag- 
netic leads by employing an equation of motion method. When the ferromagnetic leads are coupled 
with parallel spin polarization, we find three peaks in the single-particle excitation spectra. The 
middle one is the Kondo resonance caused by the orbital degrees of freedom. In magnetic fields, 
the Kondo effect vanishes. However, at a certain magnetic field new two-fold degenerate states 
arise and the Kondo effect emerges there. In contrast, when the ferromagnetic leads are coupled 
with antiparallel spin polarization, the Kondo effect caused by the spin (orbital) degrees of freedom 
survives (is suppressed) in magnetic fields. We investigate the field dependence of the conductance 
in the parallel and antiparallel spin polarizations of the leads and find that the conductance changes 
noticeably in magnetic fields. 

PACS numbers: 72.15.Qm, 73.63.Kv, 72.25.-b, 85.75.-d 



I. INTRODUCTION 

The Kondo effect, which was originally studied in di- 
lute magnetic alloys, is a typical phenomenon caused 
by local electron correlation. Recent progress on nano- 
processing techniques has enabled us to observe the 
Kondo effect in quantum dot (QD) systemsi*^. Since the 
QD system has many tunable parameters, various aspects 
of the Kondo effect can be controlled. For instance, the 
multiorbital Kondo effect realized in vertical QDs^^ and 
carbon-nanotube QDs^:^ has attracted much interest. In 
particular, in a vertical QD, the confinement potential 
yields multiply degenerate single-particle QD energy lev- 
els called the Fock-Darwin states'""^, which are regarded 
as effective orbital degrees of freedom. The degeneracy 
of the QD states in the confinement potential can be con- 
trolled with external magnetic fields. The introduction 
of the orbital splitting of the QD energy levels realizes 
a crossover from the SU(4) to the SU(2) Kondo effect, 
which results in characteristic transport propertiesi2r2i. 
In the SU(4) Kondo effect, a higher Kondo temperature 
than that in the SU(2) Kondo effect was observed in 
experiments^^—. 

Spin-dependent current in metallic leads is another 
interesting topic in relation to QD systems. Recently, 
there have been a number of experimental studies of QD 
systems coupled to ferromagnetic leads in the context 
of spintronicsi^"~. The splitting and restoration of the 
Kondo peak has been observed in experiments on a Cgo 
QD systemiS. and a semiconductor QD system, which are 
coupled to ferromagnetic leads^^'^^. There have been ex- 
tensive theoretical studies on a single-orbital QD coupled 
to two ferromagnetic leadsi^r^. When the spin polariza- 
tion of the two ferromagnetic leads is parallel, the Kondo 
effect caused by the spin degrees of freedom is suppressed 
even in the absence of magnetic fields. On the other 



hand, the Kondo effect remains when the spin polariza- 
tion of the two ferromagnetic leads is antiparallel^^"— . 
These results motivated us to investigate a multiorbital 
QD coupled to ferromagnetic leads. One can expect that 
the interplay between the spin polarization and the multi- 
orbital effects causes characteristic features of the Kondo 
effect and the related transport properties. For exam- 
ple, in a two-orbital QD coupled to ferromagnetic leads, 
the appearance of the underscreening Kondo effect was 
discussed'^!. In carbon nanotube QD, the influence of 
ferromagnetic leads on the shot noise was analyzed^^. 

In this paper, we investigate the Kondo effect of a 
two-orbital vertical QD coupled to two ferromagnetic 
leads. By employing an equation of motion (EOM) 
method2£i2^i^i2£iHi2i, we calculate the single-particle ex- 
citation spectra (SPES) of the localized electron in the 
QD and the conductance. The SPES for several spin 
polarization values are discussed in connection with the 
splitting of the QD energy levels obtained by the poor 
man's scaling approacb^SiSli^. We also investigate the 
effects of the orbital splitting caused by magnetic fields 
on the SPES. The conductance is shown as a function of 
the orbital splitting. 

This paper is organized as follows. In Sec. HIl we out- 
line the model and the methods. We consider a QD with 
two orbitals coupled to two ferromagnetic leads. The 
system is modeled as a two-orbital Anderson impurity 
Hamiltonian where the conduction electrons have a spin 
dependent density of states (DOS). We outline an EOM 
method and a poor man's scaling approach. In Sec. IIIIl 
we show the results. In the SPES for the parallel spin 
polarization of the leads, we find that three peaks ap- 
pear and the middle one is located close to the Fermi 
energy. It is considered that this is a Kondo resonance 
caused by the orbital degrees of freedom. When a mag- 
netic field is applied, the energy levels of the QD split and 
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the Kondo effect disappears. However, in a certain mag- 
netic field two energy levels cross and the Kondo effect 
emerges again. We discuss the conductance as a func- 
tion of the energy-level splitting of the QD. A summary 
is provided in Sec. |IVl 
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A. Two-orbital quantum dot coupled to 
ferromagnetic leads 

We consider a single QD with two-orbital degrees of 
freedom coupled to two ferromagnetic leads. The system 
is described by the following Hamiltonian, 

i — L,B. k,l,(T 1,(7 

where cll\^ annihilates (creates) a conduction electron 
with wave number k, spin cr, and orbital I in the lead 
i {=L, R), £/ is the QD energy level for orbital I, d'ljj 
annihilates (creates) a localized electron with spin a and 
orbital I in the QD, and ni„ = dj^di^. Intraorbital and 
interorbital Coulomb interactions are expressed by U and 
[/', respectively, and Vikia- is the tunneling amplitude be- 
tween the leads and the QD. We set U and J7' ^ as 
realized in many QD systems, so that the double occu- 
pancy of the QD electron is forbidden and the electron 
number is at most unity. Therefore, we can neglect the 
effects of the exchange interaction. We assume that the 
orbital states in the QD hybridize with the corresponding 
conduction channels in the leads. Although the assump- 
tion of multiple conduction channels is nontrivial, it is 
known that this is relevant for certain systems, e.g. ver- 
tical QD systems and carbon nanotube QD systems^"—. 
The spin polarization of the ferromagnetic leads is de- 
fined by p = (pcit - Pcii)/{Pcit + Pczi), {0 < p < 1), 
where pda- is the spin dependent DOS of the conduction 
electron in the lead i {=L, R), which is assumed to be 
independent of orbital I. We use the energy- independent 
constant pda, which is reasonable in the wide-band limit. 
For a parallel magnetic (P) configuration, the relative 
orientation of the magnetic moments in the two ferro- 
magnetic leads is parallel as shown in Fig. [IKa). Both 
DOSs in the leads L and R are defined by pcLa — PcRa = 
[1 + (Sa,^ — Scr,i)p] po with po being the DOS of the leads 
at p = 0. For an antiparallel magnetic (AP) configura- 
tion, as shown in Fig. [TJb), the DOSs in the leads L 
and R are defined by pcLa = [1 + {^c,t ~ ^o-.i)p] Po and 
PcRa = [1 - - Sa,i)p] Po, respectively. 

In a vertical QD system, the single-particle energy level 
of the QD in a magnetic field can be described as the 



FIG. 1: (Color on-line) Schematic diagram of the magnetic 
configurations in our QD system. The relative orientations of 
the magnetic moments in the two ferromagnetic leads are (a) 
parallel and (b) antiparallel. 

Fock-Darwin state. The energy levels of two orbitals are 
split into two levels, Si and 62, in magnetic fields. We 
model this splitting as 

ei^eo + iS,,i-Si^2)Krb, (^ = 1,2), (2) 

where Eq is the center of the QD energy level and Aorfc is 
the splitting width. We call Aorb the orbital splitting in 
the following. Because of the small g factor in semicon- 
ductors, the Zeeman splitting of the semiconductor QD 
system is much smaller than Aorfc in a magnetic field, so 
that we can ignore the effects of the Zeeman splitting^S. 
We thus obtain the field dependence via the Ao, b depen- 
dence using the Fock-Darwin energy diagrant^. 

B. Equation of motion method 

We calculate the retarded Green's function G[^(i) — 
~i9{t){{di„{t), <i|g.(0)}) for a localized electron of the QD. 
For this purpose, we employ an EOM method. We per- 
form the straightforward EOM calculations until the sec- 
ond iterative procedure. In the third iterative EOM 
procedure, we use the approximation which successfully 
yields the closed equations of the higher-order Green's 
functions in a single-orbital QD system^^i^. In addition 
to this approximation, we include the higher-order terms 
expressing the scattering processes of the electrons in the 
same orbital with different spins, and in the different or- 
bitals with the same and different spins. We neglect the 
higher-order terms expressing other higher-order scatter- 
ing processes. Details of the calculations are described in 
Appendix. 

In the condition of the strong interactions U and U' ^ 
\ei\, we obtain the equations for the Green's function as 

cr (, A = 1 ~ {ni^) ~ {n-J - (nj-) 

' --e.-s?,(-)-sLH-sLM-sf,H' 

(3) 

where 

E Ei^^^^'^i'^— ' (4) 
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with being the opposite spin state to a and I being the 
different orbital state from I. The Fermi energy is set at 
uj — Q. The self-energy Yj]^{uj) expresses the scattering 
process of the electrons in the same orbital with differ- 
ent spins. The self-energies E^^(aj) and Tif^{oj) express 
the scattering processes of the electrons in the different 
orbitals with the same and different spins, respectively. 
Accordingly, spin fluctuations are introduced via Ti]^{lu) 
and Yif^(uj), while orbital fluctuations are introduced via 
I]fg.(a;) and Ef^(a;). These scattering processes play es- 
sential roles in the appearance of the Kondo effect caused 
by spin and orbital fluctuations. We thus argue that the 
present EOM approximation captures the essentials of 
the Kondo effect in the two-orbital QD. 

The electron number {nia) is determined using the 
equation. 



{nia) 



Im[GL(c.)]/(c.)dw. 



(8) 



In this way, we complete the self-consistent loop which 
consists of Eqs. ([3])-®. 

To improve the numerical results quantitatively, we 
have further introduced the second self-consistent loop. 
By analogy with the self-consistent procedure used for 
a single-orbital QD coupled to ferromagnetic leads^S, we 
have replaced the bare QD energy splittings in Yi\^(uj), 
Sf^(w) and Sf^(w) by the renormalized ones Ae[^, Ae^^^ 
and Ae^^, respectively, as 



AeL 



(9) 



(10) 



Green's function in Eq. ([3]) becomes zero. This fact sug- 
gests that the additional self-consistent procedure pro- 
vides us with the important information about the renor- 
malized energy level of the electrons in the QD. For the 
single-orbital case in the P configuration, the splitting 
of the Kondo resonance was obtained using this EOM 
method^*^. This feature was confirmed by using the nu- 
merical renormalization-group calculation^^. The split- 
ting of the Kondo resonance caused by the ferromag- 
netic leads is often characterized by the so-called local 
exchange field^"*. 

Substituting the self-consistently obtained Green's 
function into the relation 



(13) 



we obtain the SPES for the localized electron of 
the QD. Note that the spectral weights obtained 
by the EOM method are rather underestimated 
and are quantitatively not accurate, although the 
spectral profiles are qualitatively correct. The 
bias-linear conductance is expressed as^i^i^ G = 

-(27reV/i) / d^/'(^) Y.ia pU^)Wal (rf. + r« ) , 

where f'{ijo) = df{uj)/duj with /(w) being the Fermi 
distribution function and F^^^' is the hybridization 
between the left (right) lead and the QD. Note that this 
expression for the conductance is valid, if the Green's 
function G'[^(a;) is diagonalized with respect to I and 
a. When we consider the QD system with off-diagonal 
contributions as well as plural electrons, we have to also 
improve the present EOM formulation. Such a situation 
is beyond the current interest. 

In the calculation, we set a flat band and as- 
sume that Vikia is independent of i, k, I, and a 
as Vikia = V. By these reasonable simpliflcations, 
we obtain Ff, = [1 + (S,^ - 5„i)p]ro/2 and Tj^ = 
[1 ± {Sa-f — Scri)p] Fo/2, where ± denotes the P(-l-) and 
AP(— ) conflgurations, and Fq — SirpoV"^ is the total hy- 
bridization of the left and right leads at p 0. The 
conductance takes the form 



G = -2tt— 
h 



[1 + ('^.t 
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Sai)p\ , P conflguration. 



AP configuration. 



(14) 



Asf. 
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£lo 



(11) 



ei, =ei+ M^U^la) + ^]a{eia) + ^lisia) + ^A^la)]. 

(12) 

where Ae^^^ denotes the difference between the QD lev- 
els of the same orbital with different spins, and Ae^^ 
and Ae^^ denote the difference between the QD levels 
of a different orbital with the same and different spins, 
respectively, u — sia determined self-consistently 
in Eq. ((HI), the real part of the denominator of the 



C. Poor man's scaling approach 

We evaluate the splitting of the renormalized QD en- 
ergy levels for the P configuration using a poor man's 
scaling approac h^^'^^ which is known as a powerful tool 
to study the renormalization effect in the Anderson 
model^. We reduce the energy scale of the half width 
of the flat band D to where charge fluctuations are 
quenched. By integrating out the excitations in the en- 
ergy range between Di and D, we obtain the renormal- 



ized QD energy levels £;^™' and e\ 



-li 



Because of this 
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renormalization, the degenerate two energy levels yield 
the spin splitting, which is given by 



Aiscai = eff' - ^T' = -Proln^ , (15) 



D 



where Di is determined from the condition that the lower 
band edge —Di coincides with the lower renormalized 
QD level eff"' according to the procedure in Refs.— i^l: 



Ztt JJl 



(16) 



We obtain Di for given p and D from Eq. (|16l) . Substi- 
tuting the obtained Di into Eq. ([T5|). we determine the 
p dependence of Aiscai ■ 



III. RESULTS AND DISCUSSIONS 

On the basis of the formulations obtained in Sec. |lTl we 
carry out the numerical calculations, keeping a vertical 
QD and Ni ferromagnetic leads in mind. We use Tq in 
units of energy. We adopt D/Tq — 100, eo/To — —2, and 
temperature T/Fq — 0.005. Since the spin polarization 
of the Ni ferromagnetic lead is p ^ 0.3^^, we direct our 
attention to the p = 0.3 system and the p dependence up 
to around p = 0.3. 
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Fig. 2: (Color on-line) Single-particle excitation spectra (a) 
p{uj) = Pt(^) + Pli^), (b) Pri^^), and (c) Pi{uj) for several p 
values in the P configuration. The parameters are Aorb/^o = 
0, T/Fo = 0.005, and eo/To = -2. Inset: Expanded scale of 
p{uj). 



A. Spin polarization dependence of Kondo effect 

We discuss the effects of the ferromagnetic leads in the 
P configuration on the Kondo effect for Aorb/^o = 0. In 
Fig. [21 we show p{uj) = Y.aP'^i^)^ Pt('^)' ^nd pi{uj) for 
several values of the spin polarization p. Note that the 
Fermi energy is uj/Tq — 0. Since we find that piaii-^) are 
independent of Z, we set Pa-{^) — J2i Piai'^)- I^i Fig- Ula), 
for p = one conspicuous peak appears at lo/Tq ~ 0.02, 
which results from the spin-orbital symmetric SU(4) 
Kondo effect^i^-ii for U and U' » |ei|. In the SU(4) 
Kondo system, the Friedel sum rule ensures that the peak 
of the Kondo resonance appears at w/Fq ~ T^^^^'^Fq, 
where jg the Kondo temperature^. From the nu- 

merical results, we evaluate T'^'^'-^-'/Fo ^ 0.02 for p ~ 
and Aorb/Fo = 0. This value is larger than the used tem- 
perature in the calculation (T/Fq — 0.005) and thus the 
SU(4) Kondo effect appears for T/Fq = 0.005. 

For p 7^ 0, the Kondo peak splits into three peaks. 
The splittings can be clearly seen for p > 0.1 as shown in 
Fig. [2lja). The energy difference ^oj/Fq between the right 
and left peaks becomes large with increasing p, while the 
middle peak remains significantly at w/Fq ~ for finite 
p. In Fig. [2Kb) and (c), we find that the left and right 
peaks in p{uj) originate from p^{uj) and Pxi'^)^ respec- 
tively. Accordingly, it is considered that for finite p in 
the P configuration the Kondo effect caused by the spin 
degrees of freedom (the spin Kondo effect) is suppressed 




Fig. 3: (Color on-line) (a) Energy splittings between up- and 
down-spin QD levels as functions of p for the P configuration. 
Aseom and Aiscai are obtained by the EOM method and the 
scaling approach, respectively, (b) Energy difference Sui/Vo 
between the peaks of p(w) in the positively and negatively 
high-energy region as a function of p for the P configuration. 
For both (a) and (b), we set Aorb/To = 0, T/Fo = 0.005, and 
eo/Fo = -2. 



as a result of the local exchange field, while the Kondo ef- 
fect caused by the orbital degeneracy (the orbital Kondo 
effect) survives. 

We investigate the origin of the splitting in terms of 
the poor man's scaling approach. In Fig. [SJa), we 
show the spin splitting of the renormalized QD energies 
{Asscai /ro) obtained by the poor man's scaling as a func- 
tion of p. We compare it with the spin splitting of the 
renormalized QD energies {Aeeoml^o) obtained by the 
EOM. According to Eqs. © and ©, Aieoml'^o can be 
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Fig. 4: (Color on-line) QD electron numbers, (rii^), (712^), 
(nil), {n2i), and (ntot) as functions of p in the P configura- 
tion, (ntot) = J2ia(^i'^)- The parameters are Aorb/To = 0, 
T/To = 0.005, and eo/To ~ —2. Inset: (ntot) as a function of 
eo/Fo with p = 0.2. 
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defined as 

Aieom = Ail ~ ^e'f (17) 

We find that both calculations are in good quantitative 
agreement, although Aiscai/^a deviates slightly from the 
linear p dependence of Aieom/^o in P 0.2. In Fig. [2Ib), 
we show the energy difference Slo/Tq between the right 
and left peaks in Fig. [2ja) as a function of p. We find that 
Slu/Tq is proportional to p and agrees with Aeeom/ro 
quantitatively. Therefore, in the P configuration spin- 
dependent charge fluctuations induce the spin splitting of 
the renormalized QD energy, which results in the linear 
p dependence of Suj/Tq shown in Fig. ^h). 

We next investigate the electron number of the QD 
in the P configuration. In Fig. 21 we show (nj^) and 
{ntot) — J2iai''^io-) functions of p. For p = 0, the 
{nia) values are the same irrespective of / and a. As 
p increases, {ni-\-) increases, while (ni^) decreases. For 
the up-spin (down-spin) electrons, the hybridization F^ 
(r|) is enhanced (suppressed) with increasing p, which 
results in an increase (decrease) in (ni^) ((n^)). On the 
other hand, we find {nia) = {n2a) for p 7^ 0. This result 
indicates that the orbital degeneracy exists even for finite 
p in the P configuration, which yields the orbital Kondo 
effect as shown in Fig. ^a). Details will be discussed 
in the next subsection. The total QD electron number 
{ntot) is almost constant irrespective of p and approaches 
unity with decreasing Sq /Tq as shown in the inset of Fig. 
3) This property originates from the strong Coulomb 
repulsions U and U' » |e;|. 

As shown in the inset of Fig. HJa), the hump structure 
around oj/Tq = — 1 caused by charge fluctuations is al- 
most independent of p. The energy at which the hump 
structure shows its maximum is shifted from the bare 
QD energy level eo/To — —2 toward the Fermi energy 
because of the correlation effects. 
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Fig. 5: (Color on-line) Single-particle excitation spectra p(t<j) 
for the several values of Aorb/Fo in the P configuration. The 
parameters are p = 0.3, T/Fq — 0.005, and Eo/Fq = —2. 



B. Orbital splitting effects 

We investigate the effects of the orbital splitting in the 
P and AP configurations. In Fig. [SJ we show the SPES 
in the P configuration for several Aorfj/Fo values at p ^ 
0.3. For Aorfc/Fo — 0, the Kondo peak caused by orbital 
fluctuations appears close to w/Fq = 0. As Ao^b/Fo is 
increased, the peak close to w/Fq = disappears, which 
indicates that the Kondo effect has vanished. Instead, 
multipeaks appear. At Aorb/^o = 0.115, we find a small 
peak at uj /Tq ~ as indicated by the arrow, showing that 
the Kondo effect is restored. When Aorb/Fo is further 
increased, the peak close to w/Fq = disappears again 
and the multipeaks appear. 

The peak positions in p{uj) are determined by the lo de- 
pendence of the self-energies T,l^{uj), T,f^{uj), and Sfg.(w), 
which reflect Ael, Aef^, and Aef^. We show iia/To as 
functions of Aorb/Fo in Fig. [HI The vertical dashed lines 
are Aorb/^o — 0.04, 0.115, and 0.16, which we have used 
in Fig. [5] 

For Aorb/Po = 0, we find eio-/Fo = £2a/To and 
£n/^o + eu/ro- According to Eqs. (1^]). (ITUl). and pT|) . 
Ae\^ and Aef^ have two different values, and Ae\^ = 0. 
Because of these properties, for Aorb/Po — 0, the Kondo 
peak caused by orbital fluctuations appears and the split- 
ting two peaks appear in p{lS). As Ao^b/Po is increased, 
the orbital degeneracy is lifted and eicr/Po (£20- /Po) in- 
crease (decrease) with the same Ao^b/Po dependence as 
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Fig. 6: (Color on-line) Renormalized QD energy levels in the 
P configuration as functions of Aorb/Fo for p = 0.3. Tiie 
vertical dasiied lines are Aorb/^o ~ 0.04, 0.115, and 0.16, 
which correspond to the values used in Fig. [S] iif/Fo and 
e24./ro cross at Aort/Fo = 0.115. 



shown in Fig. [SI Accordingly, for Aorb/^o 0; ^^'^ 
Aif^ exhibit two different values depending on a and I, 
respectively, and Aef^ exhibits four different values. Be- 
cause of these features, the Kondo effect disappears for 
Aorb/^o — 0.04 and 0.16, and eight peaks appear in p{uj) 
as shown in Fig. [5] 

For Aorb/^o = 0.115, we find a noticeable feature: 
^itZ-To Siiid e24./ro cross as shown in Fig. |6] and new 
two-fold degenerate states arise there. The scattering 
between these two degenerate states can be expressed 
using the pseudospin-flip process. This process is intro- 
duced via the self-energy E^^(w), resulting in an uncon- 
ventional Kondo effect. The small peak close to uj/Tq = 
indicated by the arrow in Fig. [5] is a Kondo resonance 
caused by this scattering process. Since the respective 
gradients of eia/^o and e2a-/^Q are the same, we find 
only one crossing point as a function of Aorh/To. Al- 
though the Aorfc/To value at the crossing point depends 
on the parameters, the result does not change qualita- 
tively: The new degeneracy arises at a certain value of 
Aorh/To, which can be evaluated from the magnitude of 
the energy splitting ejj, — e;^ at Aorb = 0, and the gradi- 
ents of £io-/ro as a function of Aorb- The similar situation 
has been observed as a singlet-triplet Kondo effect and a 
doublet-doublet Kondo effec t It is considered that 
this Kondo effect can be observed in the field dependence 
of the conductance, which will be described in the next 
subsection. 

The results for the AP configuration are shown in Fig. 
[71 For Aort/To = 0, the spectral profile is the same 
as that of the SU(4) Kondo effect shown in Fig. [DJa). 
Therefore, the effects of the ferromagnetic leads vanish 
in the AP configuration for a zero magnetic field. For 
Aorh/To = 0.01, wc find a single peak and the SU(4) 
Kondo effect remains there. As Aorb/^o is further in- 
creased, the SU(4) Kondo peak splits into three peaks. 
The peaks at either end move further apart with increas- 
ing Aorb/^o, while the middle peak remains at w/Fq ~ 
irrespective of Aort/^o- In the inset of Fig. [71 we 
show the renormalized QD energy levels as a function 
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Fig. 7: (Color on-line) Single-particle excitation spectra p(aj) 
for the several values of Aorb/^o in the AP configuration. 
The parameters are Aorb/Fo — 0, p = 0.3, T/Fq = 0.005, and 
Eo/Fo ~ —2. Inset: Renormalized QD energy levels in the 
AP configuration as functions of Aorb/Po for p = 0.3. 
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Fig. 8: (Color on-line) Conductance as functions of Aorb/Fo 
for p = 0.3 in the P and AP configurations. We set T/Fq = 
0.005. The local maximum appears at Aorb/Fo ~ 0.115 in 
the P configuration. 



of Aorb/To for p = 0.3. We find that ei^/Tq = in/To 
and £2t/ro = £24,/roi and that the two levels £ict/Fo 
and e2a/^o deviate with increasing Aorb/^o- The re- 
sults indicate that in the AP configuration the orbital 
Kondo effect is suppressed for finite Aorb/^o, while the 
spin Kondo effect remains. 



C. Conductance 

Substituting p(w) for various Aorb/Fo values into Eq. 
we obtain the Aorb/^o dependence of the conduc- 
tance G. The results for p = 0.3 in the P and AP config- 
urations are shown in Fig. [H In the P configuration, the 
conductance first decreases with increasing Aorb/^o be- 
cause of the suppression of the Kondo effect. As Aorb/^o 
increases further, the conductance increases nontrivially 
and takes a local maximum at Aorb/^o ~ 0.115. This 
nontrivial local maximum is caused by the Kondo effect 
indicated by the arrow in Fig. [51 By contrast, in the 
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AP configuration, the conductance has its maximum at 
Aorb/^o ^ 0.01. We discuss the origin of this maximum. 
For Aorb/Fo < 0.01, the SU(4) Kondo effect remains as 
shown in Fig. [71 However, the spectral weight close to 
the Fermi energy at Aorb/^o — 0.01 becomes slightly 
larger than that at Aor b/Fo = as shown in Fig. [71 For 
further increase in Aorb/Fo, the orbital Kondo effect is 
suppressed. The contributions of two spectral peaks at 
both ends to the spectral weight at uj /Fq ~ are reduced 
with increasing Aorb/^Oj leading to a monotonic decrease 
in G. Accordingly, the conductance has a maximum at 
Aorb/To - 0.01. 

We now discuss the effects of the polarization p on the 
Kondo temperature. In the AP configuration, the SU(4) 
Kondo effects at Ao,.b/Fo = is considered to appear 
even for finite p and y^^(*) is independent of p. On 
the other hand, when p is introduced for Ao^b = in 
the P configuration, the SU(4) Kondo effect is destroyed 
and the SU(2) orbital Kondo effect survives as shown in 
Fig. mja). We evaluate the Kondo temperature from the 
inflection point of the temperature dependence of G in 
the P configuration (not shown in Figure). For p — 0, 

we obtain T^^^*^ /To - 0.02 for the SU(4) Kondo ef- 
fect. This value is consistent with that obtained from 
the peak energy of the Kondo resonance shown in Fig. 
[la). For p = 0.3, we obtain T^^^'^^/Tq - 0.005 for the 
SU(2) Kondo effect. These results are reasonable, be- 
cause T^^^^ is expected to be higher than T^^'^K Ac- 
cordingly, the effective temperature scaled by T^^^^ for 

p = is lower than that scaled by Tj^ for p = 0.3 
in the P configuration. This result qualitatively explains 
the conductance at Aorb/Fo = shown in Fig. [SI be- 
cause in the AP configuration the SU(4) Kondo effect 
appears at Aorb = 0: At Aorb/Fo = 0, the effective tem- 
perature in the P configuration is higher than that in the 
AP configuration, which yields the smaller G in the P 
configuration. It is known that the EDM tends to un- 
derestimate the spectral weights and thus the magnitude 
of the conduct ance^Sii^. Actually, the conductance for 
p = and Aorb = does not approach the value in the 
unitary limit 2e^//i for the SU(4) Kondo effect. To obtain 
more quantitatively accurate results, another numerical 
calculation by, for example, a numerical renormalization- 
group may be effective. Such a research is our future 
study. 

Finally, we discuss the conductance of a vertical QD 
coupled to two Ni ferromagnetic leads. To develop a 
quantitative discussion, we use Fo = 30 meV, which 
was estimated in a Ggo experiments^. The range of 
Aorb/Fo = 0.2 in Fig. [H thus corresponds to Aorb — 6 
meV. According to the Fock-Darwin state with an elec- 
tron confinement energy of 3 meV-, we evaluate that 
Aorb = 6 meV corresponds to the magnetic field B ^ 
3.5T. On the basis of these results, we discuss the 
change in the conductance in external magnetic fields^. 
We first set the QD coupled to the Ni ferromagnetic 
leads in the AP configuration. The conductance in a 



zero magnetic field is G\b=q ^ 0.952 e^/h as shown 
in Fig. [51 When the magnetic field is increased, the 
AP configuration changes into the P configuration at 
B - O(10-S)Ti^. When the magnetic field is further 
increased, we observe a nontrivial local maximum at 
B - 1.9T and then obtain G|b=3.5 0.367 e^/h at 
B ^ 3.5T. We evaluate the ratio of the change in G as 
AG = (G|b=o - G|b=3.5t) /G|s=3.5T - 1.59. This value 
is nearly twice that in the single-orbital QD coupled to 
the Ni ferromagnetic leads, where the change ratio in G 
between the AP and P configurations is estimated to be 
AG ^ 0.81. Note that AG has the same meaning as tun- 
nel magnetoresistance. We have demonstrated that the 
large gain in AG can be controlled by an external mag- 
netic field in a two-orbital vertical QD coupled to two 
ferromagnetic leads. 



IV. SUMMARY 

We have investigated the Kondo effect of a two-orbital 
vertical QD coupled to two ferromagnetic leads using the 
EOM method. We have shown that in the P configu- 
ration the orbital Kondo effect remains, while the spin 
Kondo effect is suppressed. In magnetic fields, the or- 
bital Kondo effect also disappears, because the energy 
levels of the QD are completely split. However, at a cer- 
tain magnetic field, two of the four energy levels cross 
and the Kondo effect newly emerges there. This Kondo 
effect can be possibly observed in experiments. In the 
AP configuration the spin Kondo effect remains, while 
the orbital Kondo effect is suppressed in magnetic fields. 
We have demonstrated that the change ratio in the con- 
ductance between the AP and P configurations can be 
controlled by an external magnetic field. In a typical 
two-orbital QD coupled to Ni ferromagnetic leads, the 
change ratio of the conductance is larger than that in a 
single-orbital QD system. 
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Appendix A: Derivation of Eqs. ([3])- ([7]) 

We begin our discussion with the retarded Green's 
function G[^(i) = -i6l(t)({d/^(i), d|^(0)}). In the first 



8 



iteration of the EOM procedure, we obtain 



z|gLW = 6{t)+siGlit)+ E E^^ 

i=L,R k 

+ UUi,{t) + U'Ai,{t), 



(Al) 



where G[fc,i^(i) 



-^^{t){{cmAt),dl{0)}), 



-t9{t){{ni^{t)di^{t),dl{0)}), and Ai^{t) 



E^i^i)J{nmo.{t)dUt), dim). We perform the 
second iterations of the EOM procedure for G\^ i^{t), 
\ii„{t), and A/cr(0- For Glf.i^{t), we obtain the closed 

form as idG\^,,{t) / dt = e.kMk,iA^) + ^^^'-^[,(0- 
the latter two ones, the higher-order Green's functions 
appear as 



^^n^.w = Msit) + eiUi„{t) - E E^'fc'-^'Vw- E E^»*fc'-^^Vw+ E E^^fc'-^^'^.'-w 



i=L,i? k 



(A2) 



i-A,,(t) - E {nr,,c.)5{t)+eiKUt)- J2 J2 V,kiaY^^^j,it) - E E E ^^fci^^.L^,;. W 



+ E E E ymaYgma.M(t) + U E + E E (A3) 



r 



where is the 

the QD with spin 

^ikM^^^I 
yA 



electron number of 
a and orbital ^, and 

-m{{4kwit)dm{t)dUt),dlm), 

-ie{t){{dl{t)dut)c,ua(t),dlm), 

= -ie{t){{ni^{t)cMa{t),dl{0)}), 
-iOit){{n^a{t)ni^{t)di„it),dliO)}), 

-m{{clmc.it)dmc.mAt),dim), 

-I0{t) {{dl^ {t)dl, {t)Cikma (t) , dl (0)}) , 

= -t9{t) {t)cikia (t) , dl (0)}) , 
{t)dUt),dl{0)}). 
We perform the third iterations of the EOM procedure 
for these higher-order Green's functions which appear 
in the right-hand side of Eqs. (|X2|) and (|X3)) . In 
this calculation, we adopt the approximation which 
successfully yields the closed equation of the Green's 
functions in a single-orbital QD system^^i^. According 



yB (f\ 

ikma,la \ / 
^ikmaJai^) 



''man/3 ,la 



it) 



to this approximation, we neglect the following types of 
higher-order terms {{4kiai*)dia{t)ci'k'icr{t),dl{0)}) 
and {{dl{t)cikia{t)ci'k'ia{t),dia{Oy}) as well as 
the expectation values {clic^dia) and {dlciUa), 
and decouple {{cikw{t)'^ c^,k'ia(t)dia(t),diam}) ^ 
Siz'Skk'f{£zka){{dia{t),di„{Oy}). In addition to this 
approximation, we decouple the number operator of the 
localized electron and replace it by its expectation value. 
Furthermore, we turn our attention to the scattering 
processes of the electrons in the same orbital with 
different spins, and in the different orbitals with the 
same and different spins. We include the higher-order 
terms which express these scattering processes and 
neglect those expressing other higher-order scattering 
processes. Using these approximations, we obtain the 
following equations of the higher-order Green's functions 



(A4) 



^-^XFk,la{t) - ^ik-ffX^k,lAt) + Vi*klafi^ika)G^ai't) ~ Vikla^lait), 



(A5) 
(A6) 
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«^W™a,;<,(i) = {n,no){nia)5{t)+ \ ei + U + U' ^ (n^,) | WmaMit) - ^V^kla{nracc)X(li„{t) 



i^Yi,iAt) = I e^ka +ei-Sm-U' + U' J2 (''J^) + U^'^i-) 1 ^^''k,ia{t) + [/(e.fc^) - {nmc)] VmaG\,{t), (A9) 

'JtY^k,iA*) = e^.k.Yg^iAt) + {nmo.)V.MaG\,{t), (AlO) 

'^~77 ^ man I ait'} — {j^rna. ^manp .la 

it). (All) 



I 

We solve these coupled differential equations by the corresponds to ?7 and |£; |, we derive the equations 

Fourier transformation with respect to t. In the condition for the Green's function as shown in Eqs. (jS])-©- 
of the strong Coulomb interactions U and J7' — > oo which 
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